Strong-field approximation for harmonic generation in diatomic 

molecules 

\^ C. C. CHIRILA and M. LEIN 

O 

O ■ Max Planck Institute for Nuclear Physics, 

<N . 

Saupfercheckweg 1, 69117 Heidelberg, Germany 

O . 

O : (Dated: February 2, 2008) 

CN ■ Abstract 



The generation of high-order harmonics in diatomic molecules is investigated within the frame- 
work of the strong-field approximation. We show that the conventional saddle-point approximation 
is not suitable for large internuclear distances. An adapted saddle-point method that takes into ac- 



^ \ count the molecular structure is presented. We analyze the predictions for the harmonic- generation 

\ spectra in both the velocity and the length gauge. At large internuclear separations, we compare 
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I. INTRODUCTION 



The interaction of strong laser pulses with atoms or molecules leads to a variety of in 



teresting phenomena. Among them, the generation of high-order harmonics [1|, y, y, I^J is 
of particular interest, as it serves as a source of high-frequency coherent radiation, along 
with the generation of attosecondpulses. The possibility of their control and optimization 
promises fascinating applications 0,0]. 

For atoms, harmonic generation (HG) has been intensively studied and is well understood 

nn 

in terms of the three-step semiclassical mechanism j8|, |9|| , also known as the simple-man's 
model. In the framework of the strong-field approximation (SFA), the simple-man's model 
has its quantum- mechanical counterpart in the model proposed by Lewenstein et al. jlol |. 
The Lewenstein model for HG in atoms provides a good qualitative and quantitative agree- 
ment with ah initio calculations, which become prohibitive at high laser intensities. The 
advantage of the model is at least two- fold: First, in contrast to the integration of the 
time-dependent Schrodinger equation (TDSE), it provides a physical understanding of the 
underlying mechanisms and a certain degree of analytical description. Second, the HG 
spectrum can be obtained with much less computational effort than from the TDSE. This is 
essential for calculating macroscopic effects during the propagation of the harmonic radiation 
in the medium. 

The study of HG in molecules is at its early stages. Since molecules have additional 
degrees of freedom and more complicated symmetries, the physical phenomena are much 
richer. Experiments have shown that the ionization of molecules is greatly influenced by 
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the molecular orientation with respect to the laser polarization direction With the 

HG process depending greatly upon ionization, one has a similar effect for the emission of 



the harmonic radiation 



Ireatiy u 



. Moreover, some molecules can withstand laser intensities 



much higher than in the atomic case 



As a consequence, harmonic radiation can 



be extended to higher photon energies [lJ|. Furthermore, in the semiclassical picture the 
returning wave packet encounters the molecular core with several atomic sites whose contri- 
butions interfere. Thus, the multicenter nature of molecules can lead to interference effects 
in HG Il6l . Il7 l. Similarly, diffraction is predicted to occur in the above-threshold ionization 
(ATI) |l8l . |l9| spectra. By varying the orientation of the molecule, a certain harmonic range 



can be maximized or minimized 1£ , 
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A proposed mechanism to generate harmonics of very hi gh frequency involves probing 
a molecule stretched well beyond its equilibrium distance [2^. Note that pump-probe 
experiments have already been carried out, where such large internuclear distances have 
been created [2^. The large internuclear separation allows the electron to be detached 
from one site and accelerated by the laser field towards another site, where is captured 
and emits an energetic photon. According to the simple-man's model, the maximal kinetic 
energy before recombination can be as big as 8Up in such a process 24, 2^, where Up is the 
ponderomotive energy of the electron in the external field. 

The theoretical description of HG in molecules is largely based up to now on the nu- 
merical integration of the Schrodinger equation in one spatial dimension two 



dimensions 
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29 1 or in three dimensions 
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22I |30| . Other approaches include two- 



level models j30|, |3l|. The numerical calculations for molecules require more care than for 
atoms since one needs to study the dependence on orientation and bond length. Including 
the vibrational motion of the nuclei poses yet another difficulty. Hence, developing theoret- 
ical models can be a way to alleviate these difficulties and provide deeper understanding of 
the physical mechanisms. The first choice would be the extension of the Lewenstein model 
to the molecular case. As it turns out, due to the multicenter aspect of molecules, this 
is not a straightforward task. For example, the ground state required in the Lewenstein 
model is not known analytically for molecules. One has to use approximations, such as the 
linear combination of atomic orbitals (LCAO) approximation. It was shown that the ground 
state must have the correct symmetry in order for the model to be able to describe at least 
at a qualitative level the dependence of ionization on the molecular orientation. This can 

for 



be seen from the successful molecular- Ammosov-Delone-Krainov (MO-ADK) theory 
molecular ionization. Using the correct ground-state symmetry, the Lewenstein model for 
molecules seems to give good qualitative predictions j33[ for the alignment dependence of 
the HG yiel ds. For a more detailed discussion on the importance of the orbital symmetry, 
see Ref . 1^4 1 . 

The choice of gauge in analytical approximations is known to be difficult. While results 
from the TDSE are invariant under gauge transformations, this is not so for the SFA (for a 
discussion in the atomic case, see Ref. 35|). The SFA was used in 25 1 for a three-dimensional 
zero-range potential with two centers. Considering the harmonic emission from an atom 
which is displaced from the origin of the coordinate system by an arbitrary distance, the 
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authors show that spurious effects, hl<;e the presence of even harmonics in the spectrum arise 
in both the velocity and the length gauge. While in the velocity gauge the odd harmonics 
are invariant under the spatial translation, in the length gauge this invariance holds for 
neither the even nor the odd harmonics. Moreover, for the calculation of the harmonic 
spectra there is a choice between dipole, momentum and acceleration formulations. In 



Ref. 



22| |. a simplified version of the Lewenstein model is adapted to the molecular case. 



The authors show that the two-center interference minima can be well reproduced provided 
that the acceleration formulation is used to calculate the HG spectrum (i.e., the dipole 
acceleration is calculated by means of the Ehrenfest's theorem, as the expectation value of 
the force acting on the active electron). The dipole formulation predicts interference minima 
in severe disagreement with the TDSE results. 

In the present paper, we present a detailed analysis of the Lewenstein model for molecules, 
investigating both the length and the velocity gauge. We show that the conventional saddle- 
point approximation used for the Volkov evolution operator leads to unphysical results for 
large internuclear distances. We propose a modified version of the saddle-point approxima- 
tion, taking into account the structure of the molecule. As a consequence, one obtains the 
HG amplitude as a sum over different electron trajectories, with ionization and recombina- 
tion taking place at different molecular sites (the "exchange term" from Ref. 25|) or identical 
ones ("the return" term from Ref. ^^). We show that only in the velocity gauge, the new 
saddle-point formulation gives agreement with the cutoffs predicted by the simple-man's 
model S^- 

The paper is organized as follows. In Sec. |n] we revisit the conventional Lewenstein 
model and its extension to molecules, with the molecular ground state approximated as a 
linear combination of atomic orbitals. Both the length gauge and the velocity gauge are 
studied comparatively throughout this work. The resulting expressions are then analyzed 
within the conventional saddle-point approximation in Sec. IIIH and an adapted version of the 
saddle-point approximation is proposed. This new version takes into account the molecular 
structure, which now appears explicitly in the emerging electronic trajectories derived from 
the saddle-point equations. Sec. IIVI is devoted to the numerical evaluation of the different 
analytical expressions of the HG amplitude. The harmonic cutoffs for different internu- 
clear distances are compared to the cutoffs predicted by the simple-man's model applied 
to molecules [2^. It is shown, both analytically and numerically, that only the velocity 
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gauge together with the adapted saddle-point approximation yields cutoffs in agreement 
with the simple-man's model. Finally, Sec. |3 contains a short summary, our conclusions 
and perspectives. 



II. THEORETICAL APPROACH 

A. The strong-field approximation (length gauge) 

To obtain the strong-field approximation for the harmonic-generation process, we fol- 
low closely the approach of the Lewenstein model of HG in atoms jl^. We start from the 
dipole approximation and the length-gauge Hamiltonian for a molecule with fixed ori- 
entation irradiated by a laser field linearly polarized along the x axis (atomic units are used 
throughout): 

H = -^ + V{v,K) + E{t)x, (1) 

where V is the Coulomb interaction with the two protons and R is the internuclear distance, 
here taken as a parameter. Its orientation with respect to the laser polarization direction 
may vary. The molecule is in the (x, y) plane and the laser propagates along the z axis. 
We label the two nuclei as nucleus A at the position — R/2 and nucleus B at position R/2, 
respectively. 

Following Ref. we assume that (a) no bound states other than the ground state 
are populated, (b) the depletion of the ground state can be neglected, and (c) while in the 
continuum, the active electron does not interact with the core. Although we focus on the 
hydrogen molecular ion, for molecules such as the hydrogen molecule, one may additionally 
assume that only a single electron can become active, i.e., if one of the electrons has been 
excited into the continuum, the second electron will not couple to the field and will always 

remain in the lowest electronic state of the molecular ion. The expression for the electronic 

□ 

dipole moment is the same as given in for an atom: 
t 

D(t) = -i J dt' Jd^p d;^^[p + A{t)]dion[p + Mi'), t'] expH^(p,t',t)] + c.c, (2) 



with S = //, (it"{[p + A(t")]^/2 + Jp} being the semiclassical action and Ip the ionization 
potential of the electronic ground state. (The dependence on the internuclear distance R 
was dropped in Eq. (0) for clarity. In the following we make this dependence explicit.) Here, 



A(t) = — E{t')dt'. The HG spectrum of light polarized along a certain direction e is 
obtained by modulus squaring the Fourier transform of the dipole acceleration along that 
direction: 

e- a{n) = dte -Bit) expiint), (3) 
Jo 

where the integration is carried out over the duration of the laser pulse, Tp. Due to the 
anisotropy of the molecular system, in contrast to atoms, the emitted radiation can be 
polarized along other directions than the laser polarization axis. Here we consider only the 
harmonic radiation polarized along the direction of the laser electric field, x. 
In Eq. ((21), dion is the ionization amplitude 

(iion(k,R,t) = (^y(k)|E(t)x|^o(R)). (4) 

It has the simple interpretation of the electron transition amplitude from the ground state to 
a Volkov state |?/'y(k)) of momentum k. For the calculation of dion in © and dj-cc below, we 
use only the spatial part of the Volkov solution. The temporal phase of the Volkov solution 
is included in the semiclassical action. 
For the recombination step, we have 

d,ec(k,R) = (V^y(k)|-r|^o(R)), (5) 

whereby an electron under the influence of the external field only, described by a plane- 
wave Volkov solution, recombines with the molecular core. Because we employ plane waves 
to describe the electron motion in the continuum, some peculiarities in the predictions for 
the harmonic spectra are expected. The plane-wave description is more accurate for fast 
electrons returning to the ionic core after spending a longer time in the continuum. As 
a consequence, the predictions are expected to be more accurate for high-order harmonics 
than for low-order ones, where other effects not included here, such as the influence of the 
binding potential, come into play. A shortcoming of the current SFA formulation is that 
the ionization and recombination amplitudes are not translationally invariant, due to the 
fact that the Volkov wave function and the ground state are not orthogonal. For atoms, the 
predictions of the SFA for harmonic generation are in good agreement with the ab initio 
results, provided that the atom is considered in the origin of the coordinate system. For 
molecules, the nuclei cannot be both in the origin of the coordinate system and the lack of 
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translation invariance causes more serious problems. We discuss these in the later sections 
of this work. 

In order to calculate the dipole moment, one needs to know the ground state of the system. 
For the H-atom, this is known analytically, but for molecules one has to use approximations 
to obtain an analytical form. Here we adopt the linear combination of atomic orbitals 
(LCAO) approximation, i.e., the molecular ground state is taken to be: 

7/>o(r, R) = , ^ [M^ + R/2) + ^lJh{v - R/2)] , (6) 

with iphij) being the hydrogen ground state. Furthermore, s{R) = exp(— i?)(3 + 3i? + i?^)/3 
is the overlap integral between the two atomic orbitals, needed to ensure normalization. By 
using Eq. (jHl) in (jlj) and the ionization amplitude reads 

. /k-R\~,, , fk-R\dMky 



1 + siR) 



(7) 



where is the Fourier transform of the hydrogen ground state wave function, iphik) 

(7r\/2)^^(k^/2 + 1/2)^^. The recombination amplitude has a similar form 



drec(k,R) 



1 + s{R) 



R ^ fk-R\ ~^^^^ , _ fk-R\ dtph(k) 

dk 



R . /k- R\ ~ , 

"2" ^^"^ \~^2~ j'^''^ ' + 



(8) 



In comparison to the atomic case, the transition amplitudes ((Zj) and (jH)) depend now on the 
internuclear distance R and the molecular orientation. Thus, two-center interference effects 

Before proceeding to the next section, we briefly compare the exact bound state to the 
LCAO form that we use in the ionization and recombination amplitudes. The internuclear 
distance chosen here is the equilibrium distance in Hg". Fig.^shows the comparison between 
the gradient along the internuclear axis of the Fourier transform V'o(k) of the electronic 
ground state, obtained by using the LCAO wave function and, correspondingly, by using 
the 'exact' one. The dips in the figure appear at zeros of the two functions. The projection 
of the gradient on the internuclear axis is useful when studying molecules aligned with the 
laser field. In this case the electron momentum is mainly directed along the molecular axis 
and the recombination amplitude relates to the quantity showed in the figure. One can see 
that the LCAO approximation works well. 
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FIG. 1: The absolute value of the projection on the internuclear axis of VV'o(k) is shown for the 
LCAO approximation (continuous curve) and the exact ground state (dashed curve) of The 
internuclear distance is -R = 2 atomic units. 

B. The strong- field approximation (velocity gauge) 

To complete the analysis of the strong-field approximation, we adopt as well the velocity 
gauge. The Hamiltonian reads: 

g J-'^^^W'% ^(,R). (9) 

with A(t) the vector potential. Similar to the results in the length gauge (jll A|) . the dipole 
moment is: 

t 

D(t) = -i Jdt' Jd^p d*^^{p,t)dion{p,t') exp[-iS{p,t',t)] + c.c. (10) 


The differences with respect to the length-gauge form occur in the matrix elements describing 
the ionization and the recombination. This is due to the different form of the Volkov wave 
function and of the interaction Hamiltonian. For the ionization amplitude (iion(k, R, t) = 
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(^i/(k)| - iV ■ A{t) + A2(t)/2|V^o(R)) we obtain 

cos(^)^,(k). (11) 

The recombination amplitude has the same functional form as in Eq. (jH)), but it appears 
in Eq. (fTUI) with an argument different from Eq. The difference is that in the velocity 
gauge, the drift momentum p appears in the matrix elements, instead of the kinematical 
momentum p + A{t). 

Although the length and the velocity gauge should be equivalent, significant differences 
appear in the predicted harmonic-generation spectra. Employing the saddle-point approx- 
imation, a detailed analysis of the two gauges is given in the following. We show that the 
presence of two centers makes the usual saddle-point approximation used to sum over the 
electron momenta in Eq. (jl(J|) questionable. This is due to the oscillatory terms present in 
the Fourier transform of the ground state [e.g., the factor cos(k- R/2) in Eq. (fTT|) ]. For large 
internuclear distances, it becomes necessary to take these terms into account when applying 
the saddle-point method. 



rfion(k,R, t) 



1 + s{R) 



A(t) -k 



III. THE SADDLE-POINT APPROXIMATION AND THE CORRESPONDING 
ELECTRONIC TRAJECTORIES 

A. Length gauge 

1. The case of small internuclear distance 

The dipole expression in contains an integral over all possible intermediate electron 
momenta p. The presence of the semiclassical action ^(p,^,^') makes the integrand highly 
oscillatory and therefore the saddle-point approximation can be used. Assuming that the 
rest of the integrand is a smooth function of momentum, the saddle-point momentum is ob- 
tained from the condition that the action is stationary, i.e., VpS'(p, t, i(:')|p=p^ = 0, implying 
Ps(t',t) = — // dt" A(t") / {t — t'). The dipole moment along the laser polarization direction 
becomes: 

/ r 27r 13/2 
D,it) = -t dt' C,,Jp,(t',t) + A(t)] exp{-z5[p,(t',t),t',t]} 

J .e + i(t — t']. 







3/2 

e + i{t-t'). 

xciion[p.(t',t) + A(t'),t'] + c.c, (12) 



with e a small cutoff parameter. The term (t — t') is related to the spreading of the 
electron wave packet during the motion in the continuum. 
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For a complete description in terms of complex electron trajectories j36l . 
contribute to the HG spectrum, one needs to apply the saddle-point approximation for 
both the integration over the time t' in (fT^ and over t in (jS)). This gives two more saddle 
equations, one for the tunnelling time t'^ and one for the harmonic emission time t^. As 
all the saddle equations are obtained from the semiclassical action, which is the same as in 
the atomic case, the electron trajectories do not reflect the structure of the molecule. The 
only differences lie in the expressions for the recombination and ionization amplitudes, but 
the saddle momentum which governs the trajectories in the simple-man's model does not 
'feel' the molecular structure. Although they are the same as in we repeat them here 
for completeness: 

[p.(t;,g + A(t-)]^ _ 

MU)±AM! = ^_,^. (13) 

They simply state that the electron is born in the continuum with approximately zero kinetic 
energy and that the energy is conserved when the electron recombines into the ground state, 
followed by emission of a photon with energy VL. Equations (fTT^ can be used to determine 
the cutoff in the harmonic spectrum, which is predicted to occur for a monochromatic laser 
field close to ^2 = 3.17f/p + Jp, the same as in the atomic simple-man's model. 

Unlike the H-atom case, for molecules both the ionization (|7j) and the recombination (jHJ 
amplitude contain an oscillatory factor, whose argument is proportional to the internuclear 
distance R. This means that for increasing i?, the integrand becomes more and more os- 
cillatory as a function of t and t', i.e., the above saddle-point equations are not valid and 
the cutoff of the harmonic spectrum increases. We find this effect present in the numeri- 
cal calculations. Moreover, the extra oscillatory factor makes the saddle-point method for 
momentum questionable. Hence, it becomes necessary to include the extra oscillations by 
adding them to the semiclassical action, when the saddle-point approximation is applied. 
We discuss this approach and its consequences in the following. 
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2. The case of large internuclear distance 



To account for the extra oscillatory factors, we include them in the oscillatory function 
used to derive the saddle momentum. In the simpler method of the previous section, only the 
oscillations due to the semiclassical action were taken into account. The present procedure 
reveals a more accurate picture of the harmonic generation process. It changes the definition 
of the complex electron trajectories, which now depend on the structure of the molecule. We 
proceed by re-writing the ionization and recombination amplitudes in the following form: 



From the product (i*pc^x[p + A(i(:), R](iion[p + A(t'), R, t'] of the two amplitudes (fTijl there 
are four terms arising for different combination of signs in the exponential. Each of these 
terms gives rise to a different oscillatory behavior of the integrand. Therefore, the saddle 
equation for the electron momentum is different for each of the four cases. We treat them 
separately and emphasize the physical picture they convey. 

1. The combination with the prefactor r*j^i^ leads to an oscillatory behavior through 
the phase S'(p,t',t) + p • R + [A{t) + A(t')] • R/2. The condition for the saddle 
momentum reads VpS'(p, t', t)|p=p^ = — R, i.e., the saddle momentum is Ps(t',t) = 
- dt"A{t")/{t - t') - R/(t - t'). The absolute value of the integral over A{t") can 
be estimated to be smaller than 2ao, where = Eq/uj"^ is the classical amplitude of 
oscillation of the electron in the external field. Thus it becomes apparent that the 
present modification of the saddle-point method is certainly required for R > 2ao, 
giving us a more precise definition of large internuclear distances. 

Applying the saddle condition to perform the integration over t' in Eq. and over t 
in Eq. (jS)), one obtains for the tunnelling time t'^ of the electron in the continuum 



(iion(k, R, t) = i+ exp(zk ■ R/2) + z_ exp(— zk • R/2) 



(k,R) 



r+ exp(ik • R/2) + r_ exp(-ik • R/2). 



(14) 



Ps{t',,t) + A{t')]' 



+ Jp + E(t;) ■ R/2 = 0, 



(15) 



2 



and for the emission time tg of a harmonic with frequency ^2 



[p,(t;,t,) + A(g]2 



n-Ip + E{ts) ■ R/2. 



(16) 
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The expression of the saddle momentum along with the condition for the tunnelling 
time ()15|) define an electron trajectory such that, after the electron was 'born' in 
the continuum at nucleus i? at a time close to Re(tg), it oscillates in the field and 
recombines with the molecular core at nucleus A. We call the harmonics generated by 
tunneling and recombination at different molecular sites transfer harmonics. 

Equations (|15p and are conditions for the conservation of energy during ioniza- 
tion and recombination, respectively. In these equations, the ionization potential is 
unchanged (remember that the SEA considers the ground state not dressed by the 
external field). Only the continuum is shifted by the potential energy of the electron 
in the external electric field at the two sites, A and B. This way, there is an additional 
unphysical energy difference between the Volkov state and the ground state of the 
electron, which can be seen in Eqns. (jl5|) and 

Similarly, the term with the prefactor r!ii+ describes a process where the electron is 
'born' at the nucleus A and recombines at the nucleus B. The equations describing 
this process are as given above, with R changing sign. 

2. The combination with the prefactor r'^i^ leads to the phase ^(p, t', t) + [A(t) — A{t')] ■ 
R/2. The condition for the saddle momentum reads VpS'(p, t', )f:)|p=p^ = 0, i.e., the 
saddle momentum is Ps(t',t) = — J^, dt" A(t") / {t — t'). Eor the tunnelling time one 
obtains: 

IP-"'-"+^"''l^ + - m) ■ R/2 = 0. (17) 

The expression of the saddle momentum along with (|T7jl define an electron trajectory 
such that, after the electron was 'born' in the continuum at nucleus A at a time close 
to Re(t'^), it oscillates in the field and recombines with the molecular core at the same 
nucleus (we call these direct harmonics). The saddle condition for the emission time 
ts reads 

[P3(t;,t,) + A(t,)]^ = - Jp + E(t,) ■ R/2. (18) 

Eor the ionization and recombination of the electron at the nucleus A, the ionization 
potential is unaffected by the field, but the continuum is shifted by the potential energy 
E • R/2 of the electron in the external electric field. 

The term with the prefactor r* i_ corresponds to an electron 'born' at nucleus B and 
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recombining with the same nucleus. The relations describing this process are as given 
above, with R changing sign. 

The saddle equations for the emission time, Eqns. (fT^ and (fTHj). state that the relation 
between the energy of the returning electron and the harmonic frequency of the emitted 
radiation depends on the continuum energy shifts E(t'^) • R/2 and E(ts) ■ R/2. For large 
values of R, this gives rise to unphysically high harmonic cutoffs. One can make a rough 
estimate of the cutoff for large R, obtaining of the order 0{EqR), where Eq is the electric 
field amplitude. 

In conclusion, the harmonic cutoff still increases monotonically with increasing internu- 
clear distance. In the velocity gauge, the new formulation of the saddle-point approximation 
for momenta eliminates this artifact. 

B. Velocity gauge 

1. The case of small internuclear distance 

Using the velocity-gauge expressions for the ionization and recombination amplitudes 
from Sec. IIIB| and employing the conventional saddle-point approximation for the electron 
momenta, one obtains the dipole moment along the laser polarization direction 



with e a small cutoff parameter. The saddle momentum is the same as in the length gauge. 
Again, the ionization and recombination amplitudes contain oscillatory factors that give 
rise to harmonic cutoffs that increase with increasing internuclear distance. The situation is 
completely different if we take the extra oscillatory factors into account when we apply the 
saddle-point approximation. 

2. The case of large internuclear distance 

We proceed similarly to the length-gauge calculation, re-writing the ionization and re- 
combination amplitudes in the form of Eqns. From the product d*g^^^{k, R)(iion(k, R, t) 








[Psif, t),t'] + C.C., 



(19) 
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of the two amplitudes there are four terms arising for different combination of signs in the 
exponential. Each of these terms gives rise to a different oscillatory behavior of the inte- 
grand. Therefore, the saddle equation for the momentum of the electron is different for each 
of the four cases. We treat them separately and analyze the corresponding physical picture 
in terms of electronic trajectories. 

1. The combination with the prefactor r*^i_ leads to the phase S{p,t',t) + p ■ R. The 
condition for the saddle momentum reads VpS'(p, t', t)|p=p^ = — R, i.e., the saddle 
momentum is Ps(t',t) = — /// dt" A{t") / {t — t') — R/(t — t'). Using saddle-point ap- 
proximations to perform the integration over t' in Eq. (j21) and over t in Eq. (jS)), one 
obtains for the tunnelling time t'^ of the electron 

IP-"--'');^"-'l%/, = 0. (20) 

and for the emission time tg of a harmonic with frequency Vt 

[P.(i;.i.)+A(i.)]^ = Si - /, (21) 

The expression of the saddle momentum along with the condition for the tunnelling 
time ()2()|1 define an electron trajectory such that, after the electron was 'born' in 
the continuum at nucleus i? at a time close to Re{t'g), it oscillates in the field and 
recombines with the molecular core at nucleus A. 

Similarly, the term with the prefactor rli_|_ describes a process where the electron 
is 'born' at nucleus A and recombines at nucleus B. The equations describing this 
process are as given above, with R changing sign. We note that unlike for the case 
of the length gauge, the saddle equations that describe the complex trajectories agree 
with the classical equations of motions of the simple-man's model if one sets Jp = 0, 
as is usually done when casting the complex trajectories in the language of the simple- 
man's model. The harmonic cutoff in this case can be as much as 8f/p + Jp for a 



monochromatic field 



m 



25|. This will be addressed in more detail in the next section. 



2. The combination with the prefactor r'^i+ leaves the semiclassical action unchanged. 
The condition for the saddle momentum reads VpS'(p, t', t)|p=p^ = 0, i.e., the sad- 
dle momentum is Ps(t',t) = — J^, dt"A(t")/(t — t'), like in the atomic case. For the 
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tunnelling time one obtains: 



+ 1, 



0. 



(22) 



2 



The expression of the saddle momentum along with (j^^ define an electron trajectory 
such that, after the electron was 'born' in the continuum at a certain nucleus at a time 
close to Re{t'g), it oscillates in the field and recombines with the molecular core at the 
same nucleus. The saddle condition for the harmonic emission time tg reads 



The term with the prefactor r^z_ corresponds to the same kind of trajectory, but 
for the other nucleus. The saddle equations are identical to the ones that describe 
the harmonic generation in the atomic case, and therefore have the usual cutoff of 
3.17?7p + Jp for stationary laser fields. 

In the velocity gauge, the new formulation of the saddle-point for momenta eliminates the 
artifact that the harmonic cutoff increases indefinitely with increasing internuclear distance. 

IV. RESULTS AND DISCUSSION 

For the calculation of harmonic spectra, we choose a laser electric field with two optical 
cycles turn-on and turn-off and a plateau of three constant-amplitude optical cycles. Unless 
stated otherwise, we consider a molecular ion aligned with the field, irradiated by a laser 
with a wavelength of 800 nm and intensity equal to 2 x 10^^ W/cm^. The energy of the elec- 
tronic ground state is given by Ip{R) = Vbo{R) — where Vqq is the Born-Oppenheimer 
potential for the electronic ground state (including the proton-proton interaction energy). 

We focus our analysis on the variation of the cutoff energy with the internuclear distance 
R j^,!^]- particular, we consider the cutoff related to the process in which the electron 
leaves the nucleus A and recombines at the nucleus B. As discussed in ^25^, for R = (2n -|- 
l)7rQ;o (n = 0, 1, 2 . . .) the maximal kinetic energy at B can be as much as 8f/p for a linearly 
polarized monochromatic field, where ao is large on the atomic scale (for the present laser 
parameters, Oq = 23.3 a.u.). For the maximal kinetic energy, the value of the internuclear 
distance is such that an electron detached at zero electric field at A is accelerated for half an 
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FIG. 2: The maximal kinetic energy upon arriving at nucleus B of an electron detached with zero 
initial kinetic energy from the nucleus A, as a function of the internuclear distance. Full curve, 
result for a monochromatic laser field, see Ref. Q. Dashed curve, result for a trapezoidal envelope 
(see text). 



optical cycle before reaching B. This way, the energy at recombination is much larger than 
3.17t/p, as shown in Fig. |21 

In the following, we analyze the harmonic spectra for both the length gauge and the 
velocity gauge, using the conventional and the adapted saddle-point method for integration 
over the electron momenta. 

A. Length gauge results 



To better understand the harmonic spectra as they emerge from the SFA, panel (a) of Fig. 
El shows the harmonic spectrum for the atomic case. Using the conventional saddle-point 
method, the contribution of the recombination amplitude to the dipole moment simplifies 
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to: 
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(24) 



where Rx is the projection of the internuclear distance on the polarization direction and fcret 



is the momentum of the returning electron, fcret = y 2(n — Jp). From Eq. (j^^ . it is seen that 
the behavior of the harmonic efficiency for harmonic orders bigger than the molecular binding 
potential is oscillatory, while for smaller harmonics is exponential. Thus, the predictions 
of the SFA in the energy region around the binding potential Jp are distorted. This energy 
region is not described well by the SFA model, which does not fully include the effect of the 
Coulomb potential. Another remark in view of the above is that the two-center interference 
effect in HG is seen to stem mainly from the explicit form of the recombination amplitude. 
Therefore, for a good description of the interference, one has to provide a realistic expression 
for the recombination amplitude. 

Fig. ini shows harmonic spectra for different internuclear distances. The first thing to 
notice is that both the harmonic cutoff and the harmonic intensity increase with increasing 
internuclear separation. This is due to the presence of the oscillatory factors depending on R 
in the ionization and recombination amplitudes (see the comments at the end of Sec. lIII AT| . 
The arrows in Fig. IHl show the cutoff of the transfer harmonics as predicted by the simple- 
man's model (Fig. EI). Since the present SFA formulation does not include rescattering of 
the electron from the molecular sites, we expect that the maximal harmonic order comes 
from the transfer mechanism. As seen from the figure, for large values of R, there is a 
significant difference between the SFA predictions in in the length gauge and the simple- 
man's predictions. The following figures will again display the cutoffs from the simple-man's 
model for comparison. 

The contribution from the ionization amplitude in the case when R is small depends on 
terms such as sinh(i?j. ^2Ip/2) and cosh(i?a; ^ 21^/2) , as it can be shown by expanding the 
relevant quantities in Taylor series around the ionization times from the simple-man's model. 
This explains the increase of the harmonic intensity with increasing internuclear distance. 
Thus, the unphysical increase in the harmonic efficiency is exponential in R. 

So far, we have discussed the results within the conventional saddle-point method for 
summation over the electron momenta. We turn now to the second saddle-point formulation 
for electron momenta (presented in Sec. lIII A^ . Fig.|3]plots the predicted harmonic spectra. 
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FIG. 3: Calculated harmonic spectra using the length gauge and the conventional saddle-point 
method (see text), for various internuclear distances in H^: (a) i? = (the atomic case, Ip = 1); 
(b) R = O.Svrao; (c) R = ttuq and (d) R = l.Svrao- The laser and molecule parameters are 
described in the beginning of Sec. IIVI and the arrows show the cutoffs of the transfer harmonics 
calculated from the simple-man's model. 

To assess the importance of the direct harmonics, we plot them in panel (d) of Fig. |3] Their 
cutoff is smaller in this case than the cutoff of the transfer harmonics. The artifacts regarding 
the cutoff and the harmonic efficiency are present. 

B. Velocity gauge results 

For the velocity gauge and using the conventional saddle-point method, the harmonic 
spectra (Fig. Ej) are similarly unphysical as those obtained in the length gauge. The shape 
of the spectra is slightly changed together with the values of the cutoffs. The increase of the 
cutoff with increasing internuclear distance R is present. A striking difference with respect to 
the length gauge, however, arises when one uses the second formulation of the saddle-point 
method, see Fig. (jU)). In this case, the unphysical increase of the cutoff with R is eliminated. 
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FIG. 4: Harmonic spectra (same parameters as in Fig. but calculated by using the adapted 
saddle-point method and the length gauge. The lower curve in panel (d) shows the direct harmonics 
only. The arrows show the cutoffs of the transfer harmonics calculated from the simple-man's model 
(see Fig.Hl). 

C. Comparison of the harmonic generation cutoffs with the predictions of the 
simple-man's model for large internuclear distances 

We compare the cutoffs for the transfer harmonics as predicted by the SPA using the 
adapted saddle-point method in both the length and velocity gauge with the predictions 
of the simple-man's model. Fig. [7| shows the numerical results for very large internuclear 
separations. The arrows point to the semiclassical cutoffs from the simple-man's model. 
One can see that the length gauge overestimates by far the cutoffs as well as the harmonic 
intensity. In contrast, the velocity gauge together with the adapted saddle-point formulation 
succeeds in reproducing the simple-man's cutoffs, as well as keeping a bounded magnitude 
for the harmonic intensities. 
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FIG. 5: Calculated harmonic spectra using the velocity gauge and the conventional saddle-point 
method (see text), for various internuclear distances in H^: (a) i? = (the atomic case, /p = 1); (b) 
R = O.Syrao; (c) R = ttoq and (d) R = l.Svrao- The laser and molecule parameters are described 
in the beginning of Sec. IIVI The arrows show the cutoffs of the transfer harmonics calculated from 
the simple-man's model. 

V. CONCLUSIONS AND PERSPECTIVES 



The SFA is a useful tool for studying processes in intense laser fields. It is of big advantage 
for molecules, where a solution of the Schrodinger equation can be more prohibitive than 
for an atom. There are still open questions about how to adapt the SFA for molecules such 
that its level of prediction would become as good as for the atomic case. 

Here, we have analyzed the strong-field approximation for the calculation of the harmonic- 
generation spectra in ultrashort laser pulses. Using the conventional saddle-point method 
to approximate the sum over the electron momenta in the Volkov propagator, we have 
shown that both the length gauge and the velocity gauge predict harmonic intensities and 
cutoffs that increase unphysically with increasing internuclear separation. Also, the emerging 
trajectories for the electronic complex orbits jSf 
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do not reflect the structure of 
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FIG. 6: Harmonic spectra (same parameters as in Fig. but calculated by using the adapted 
saddle-point method and the velocity gauge. The lower curve in panel (d) shows the direct harmon- 
ics only. The arrows show the cutoffs of the transfer harmonics calculated from the simple-man's 
model. 

the molecule and are the same as for the atomic case. In order to account for the presence 
of the two centers, we have introduced an adapted saddle-point method that takes into 
account the molecular structure. As a consequence, we obtain explicitly different types of 
electronic trajectories and corresponding harmonics, namely direct harmonics (generated by 
tunneling and recombination of the active electron at the same site) and transfer harmonics 
(generated by tunneling followed by recombination at a different site). The interesting 
observation is that the length gauge still predicts in the adapted saddle-point formulation 
unbounded cutoffs and harmonic intensities as the internuclear separation is increased. Thus, 
the length-gauge SEA, in its present form, appears to be inappropriate for the description 
of high-harmonic generation at large internuclear distances. In contrast, the velocity gauge 
predicts bounded harmonic intensities and cutoffs of the transfer harmonics. Only in the 
velocity gauge the cutoffs agree with the values obtained from the simple-man's model. 
Not addressed in this work is the question about the role played by the excited molecular 
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FIG. 7: Transfer harmonics (same parameters as in Fig. calculated using the adapted saddle- 
point method in the length gauge (left column) and the velocity gauge (right column). The arrows 
show the cutoffs of the transfer-harmonics spectrum, predicted by the simple-man's model. 

electronic states in the harmonic emission process and whether or not it is necessary to 
include them in the SFA. Furthermore, we have not investigated the question which gauge 
is suitable for small internuclear separations, but there are indications that the length gauge 
will be preferable in this case j^]. We believe that this work is a step towards building an 
accurate SFA for harmonic generation in molecules. 
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